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A Coupled Nonlinear Spacecraft Attitude Controller/Observer With an Unknown 

Constant Gyro Bias 

J. Deutschmann and R.M. Sanner 


Abstract - A nonlinear control scheme for attitude control of a spacecraft is combined 
with a nonlinear gyro bias observer for the case of constant gyro bias. The closed loop 
system is proven to be globally stable, with zero tracking error, thus proving a separation 
principle for the given system. The nonlinear observer incorporates persistency of 
excitation, resulting in exponential convergence of the gyro bias error. 


I. Introduction 

Combined observer-controller designs for nonlinear systems are a subject of active 
research [1,2,3]. Successful design of such architectures is complicated by the fact that 
there is, in general, no separation principle for nonlinear systems. In contrast to linear 
systems, “certainty equivalence” substitution of the states from even an exponentially 
converging observer into a nominally stabilizing state feedback control law does not 
necessarily guarantee stable closed-loop operation for the coupled systems [2,4,5]. 

In this paper we consider a restricted version of this problem, in particular the task of 
forcing the attitude of a rigid vehicle to asymptotically track a (time-varying) reference 
attitude using feedback from sensors with persistent nonzero bias errors. Specifically, 
we propose utilizing an angular velocity observer from [2,3] in a certainty equivalence 
fashion with the nonlinear control law proposed in [6], and demonstrate in this case that 
the resulting system results in stable closed-loop operation, with asymptotically perfect 
tracking. 

The proof proceeds in two steps. First, we extend the analysis of the observer in [2,3] to 
the case of constant gyro bias, and use a persistency of excitation argument to 
demonstrate that the bias estimates provided by this observer are in fact exponentially 
convergent to the true bias values. The proof in [2] uses a Lyapunov argument to obtain 
a similar result, but under the assumption that the biases themselves are exponentially 
decaying; this restriction is removed here. Second, we consider the certainty equivalence 
use of these observer estimates in the nonlinear feedback control algorithm proposed in 
[6] and show that the perturbation introduced by this strategy into the closed-loop 
dynamics can be represented as a bounded function of the vehicle states multiplying the 
observer transients. Quantifying the impact these perturbations have on the Lyapunov 
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analysis given in [6], we demonstrate that the stability and convergence properties of the 
controller are, in fact, maintained in the face of the perturbations, completing our analysis 
of the proposed methodology. 

The paper is organized as follows. Section II contains definitions of the terms used in the 
controller and observer. In Section III the nonlinear observer for the constant gyro bias is 
developed and the exponential convergence is proven. Section IV presents the nonlinear 
controller design and the proof of stability of the closed loop system and the 
convergence of the tracking errors. Section V presents simulation results, followed by 
conclusions in Section VI, 


II. Definitions 

The attitude of a spacecraft can be represented by a four component quaternion, 
consisting of a rotation angle and unit rotation vector, known as the Euler axis 


£ 
ft. 

where <f> is the rotation angle, e is the Euler axis, £ and q are the vector and scalar 

portions of the quaternion, respectively. Note that ||q|| = 1 by definition. The quaternion 

represents the rotation from an inertial coordinate system to the spacecraft body 
coordinate system. A rotation matrix can be computed from the quaternion as [6] 

R(q) = (ft 2 - £ T e)I + 2z £ T - 2qS(fi) 
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where S(e) is a cross product matrix formed from the vector s. 
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A desired target attitude is represented by the quaternion, q^ =[sj,ft d ]. The attitude 

error used in the controller is defined as a rotation from the desired body frame to the 
actual body frame and is computed according to [7] 
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Similarly, in the observer, the attitude error is defined as the rotation from the estimated 
body frame to the actual body frame as 
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where q represents the attitude state of the observer. Note that ? c = 0, rj c = ±1 indicates 
that the spacecraft is aligned with the desired attitude and similarly, e 0 = 0, rj Q =±1 
indicates that the attitude estimate is aligned with the actual attitude. 


The kinematics equation for the quaternion is given as 
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where co is the spacecraft angular velocity. The angular velocity is typically measured by 
a gyro, which can be corrupted with both systematic and random errors. In this work we 
consider only the case of systematic errors. In the case of a gyro bias, the angular 
velocity can be written as 


co = co g + b 

where co g is the angular velocity from the gyro and b is the gyro bias. An estimate of the 
angular velocity is given as <b = co g + b. The bias error is defined as the difference 
between the true and estimated bias 


b = b - b (2) 

Finally, a measure of the discrepancy between the actual and desired angular velocity in 
the controller is computed as [8] 

65 c = co - R(q c )co d (3) 

which is defined such that q c = |Q(q c )o> c . 

III. Nonlinear Observer for Constant Gyro Bias 

Following the development of [2] a state observer for the bias can be defined as 

q = ^QR T (q 0 )(ci> g + b + ke 0 sgn(n o )) (4) 

b =|£ 0 s g n (fto) 
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The gain, k, is chosen as a positive constant. The R T (q 0 ) resolves the angular velocity 
terms in the observer frame. 

Computing the derivatives of q 0 in (1) and b in (2) results in the following differential 
error equations 
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b = £ 0 sgn(rj 0 ) 


The equilibrium states are (0,0,0, ± 1,0, 0,0). Again, following the development of [2], a 
Lyapunov function is chosen as 


V = |b T b + 1- 
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Taking the derivative of V 0 , and noting that £ 0 e 0 + rj 0 rj 0 = 0, yields 


Mo q 0 2: 0 

Ho Oo < 0 


Substituting rj 0 and b from above results in 

M-tMo^o 

This establishes that b, £ 0 ,andrj 0 are globally uniformly bounded; a further application 
of Barbalat’s lemma [4] shows that, in fact, e 0 — » 0 . 

Since all the signals in the observer have been demonstrated to be bounded, ? 0 and b 
can be analyzed as 


x(t) = A(t)x(t) 


and the derivative of the Lyapunov function written as 

V 0 = -x T C T Cx < 0 


where 
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By a standard argument (e.g. Theorem 4.5 of [4]), the equilibrium point x = 0 of this 
system is exponentially stable if the pair (A(t),C) is uniformly completely observable 
(UCO). Since observability properties are unchanged under output feedback, this will be 
true if and only if the pair (A(t) - K(t)C,C) is uniformly observable, with K(t) piecewise 
continuous and bounded. With K(t) chosen as 


K(t) = 
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The matrix (A(t) - K(t)C) becomes 


A(t)-K(t)C 


0 l(rj 0 I + S(€ 0 )) 
0 0 
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The state transition matrix for (6) is 


0(T,t) 


I iI(T,t) 
0 I 


with |I(T,t) = iXQi(‘?o( a )) dcr = y| t T (q 0 (a)I + S(e 0 (a)))da . The observability Gramian 
for the pair (A(t) - K(t)C,C) is given by 


W(t + T) = j;* T O t (i, t)C T CO(x, t)dt = f T 
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and the system is UCO if there exists a T>0 and positive constants ki>0, k 2 >0 such that, 
for all t > t 0 


k, I > W(t + T) > k,I 

Using Lemma 13.4 of [5], this property is assured for (7) if Q ( (t) and Q,(t)are bounded, 
and there are positive constants T : , a, and a 2 such that, for all t > t„ 


a :l - 1,'^ Qi( T )Qi T ( T )dt > a,I 
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The upper bound is satisfied since Qi(t) is bounded. To determine if the lower bound is 
satisfied, rewrite the matrix Q|(t)Q, T (t) as 

Q,(0Q,'(t) = (?j 0 (t)I + S(E 0 (t))(n o (t)I -S(£ o (0) = I - E„(t)s 0 (t) T (9) 

Since e 0 — » 0 asymptotically, for any 5>0, there exists a Ti(6)>0 such that ||e 0 || < 5 for all 
t > t 0 +T, . Taking any 5<1 and T 2 >Ti, (9) implies that 

z T [r T2 Q.(t)Q, T (T)dT]z = (1-5)(T 2 - T, )||z|| 2 

for any z in R 6 and thus (8) is satisfied. This demonstrates the required UCO property for 
the observer, and hence establishes that e 0 and b approach zero exponentially fast. 

IV. Nonlinear Controller Design 

The complete attitude dynamics for a rigid spacecraft are given as 

Hob - S(Hco)co = u 

q = iQ(q)“ 

where H is a constant, symmetric inertia matrix and u is the applied external torque, for 
example, from attached rocket thrusters. The goal of the controller is for the actual 
attitude q(t) to asymptotically track a (generally) time-varying desired attitude q d (t) and 
angular velocity w d (t), related for consistency by q d = 4-Q(q d )to d . It is assumed that 

Wd(t) is bounded and differentiable with 6 d (t) also bounded. 

The passivity based controller of [6] utilizes the composite error metric 

s = <n c + X e c = to - to r (10) 

Where from (3), co r = R(q c )to d - ae c . Taking the derivative of (10) and multiplying by 
H results in 


Hs = Hcb - Hcb r = u + S(Hto)to - Ha r 
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where 


a r = cb r = R(q c )cb d -S(S c )R(q c )to d ->.Qj(q c )« c 


and Q, (q c ) = rj c I + S(£ c ) as defined above. With these definitions, the control law 

u = -K d s + Ha r -S(Ha))o) r (12) 

for any symmetric, positive definite matrix K D as shown in [6] produces the desired 
stability and asymptotic tracking properties. 

In the current application, the control law (12) cannot be implemented because exact 
measurements of the angular velocity © are not available. Instead a certainty equivalence 
approach is employed using the estimates d> from above, resulting in 

u=-K D s + Ha r -S(Hd>)to r (13) 

where s = e>-co r , <3 C = co - R(q c )co d , and 

a r = R(q c )cb d + S(R(q c )to d )co c -XQ,(q c )to c . 

Substituting (13) into (11), along with (10), and noting that ? = s-s = b, 
a r = [S(R(q c )co d )-XQj(q c )]b, andco c - co c = © - d> = b produces the closed-loop 
dynamics 

Hs - S(Hco)s + K d s = [-S(co r )H - HS(R(q c )to d ) + XHQ, (q c ) + K D ]b (14) 

The terms on the right hand side of (14) can be rewritten as A(q c ,co d )b . Since ||q c || = 1 
by definition and |to d || < co by assumption 

7 s sup sup|A(q e ,o d (t)l<oc 

Using the Lyapunov function V c =yS T Hs, the derivative of V c along closed-loop 
trajectories of (14) satisfies the inequality 


V. = -s t K d s + s T Ab < -k D 


+ m 


where k D is the smallest eigenvalue of Kq. Using Young’s inequality [5] on the last term 
above, V c is rewritten as 




Thus, recalling from the observer analysis that |b|| is bounded, s is also seen to be 

uniformly bounded. Similarly s is uniformly bounded, since all the terms in (14) have 
been demonstrated to be bounded. Moreover, s € L 2 since for any t > tQ 
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The integral on the right hand side of (15) is finite for any t > t 0 since ||b[| is converging 

exponentially to zero. Thus seL. x r>L 2 , s e L x , hence by Barbalat’s lemma [5] 
lim s(t) = 0. 

t— 

This establishes convergence of the composite metric s to zero; now the convergence of 
the actual attitude error s c is examined. From (10), ||s[| 2 = |jfi c || 2 + 2X65je c + X : ||e c ||\ 

Since ||e|| c is bounded by definition, seL* implies that to c e L x . This in turn 
demonstrates that \ e L x since q c =yQ(q c )5 c and both terms on the right are 
bounded. Finally, e c is also in L 2 since for any t > to 
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noting that <oJe£=2rj c , the terms on the right of (16) are finite. Thus, 
e n L 2 , e c e L x and Barbalat again implies that lim e c =0. (Note that an 

t-»x 

analogous argument establishes that limS c =0.) 


V. Simulation Results 

The spacecraft attitude controller/observer design is tested with a Matlab simulation. 
The inertia matrix is a diagonal matrix with principal moments of inertia of [90, 100, 70] T 
kg-m 2 . The desired angular velocity is ©/ = [0, 0.11, 0] rad/sec and qd T =[0, 0, 0, 1] T . 
Table I lists the initial conditions for the observer and controller, as well as the true gyro 
bias and true initial angular velocity. The gains are chosen as k=l, K D =k D I with k D =6. 

Table I. Simulation Initial values 


Variable 

Initial Value 

q 

[0, 1,0, 0] T 

q 

[0, 0, 1,0] T 

to - true 

[-. 1 ,.2,-.4] t rad/sec 

b - true 

[.05, -0.05, 0.033 ] T rad. sec 



b [0, 0, 0] T rad/sec 

Figures l and 2 show that e 0 and b , respectively, converge to zero. Figures 3 and 4 
show that e c and 65 , respectively, also converge to zero. 




Figure 2. Observer Error - Ib^ 



Figure 3. Controller Error - j|e c j|. 





Figure 4. Controller Error - 65 


VI. Conclusions 

A nonlinear controller/observer is presented for spacecraft attitude applications, given a 
constant gyro bias. The nonlinear observer is globally convergent, with identification of 
the gyro bias proven through a persistency of excitation argument. The nonlinear 
controller is a passivity based controller. The control input requires the use of the gyro 
bias estimate from the nonlinear observer. The closed loop stability properties of this 
nonlinear controller coupled with the nonlinear observer are analyzed and the system is 
found to be globally stable, leading to a separation principle for the nonlinear system. 
Given a desired target attitude and target angular velocity, this system guarantees zero 
tracking error. 


References 

[1] S. Salcudean, “A Globally Convergent Angular Velocity Observer for Rigid Body 
Motion”, IEEE Transactions on Automatic Control, Vol. 36, No. 12, December, 
1991. 

[2] H. Nijmeijer and T.I. Fossen, Lecture Notes in Control and Information Sciences, 
New Directions in Nonlinear Observer Design, Springer-Verlag, 1999, pp. 135- 
159. 

[3] R. Alonso, J.L. Crassidis, and J.L. Junkins, “Vision-Based Relative Navigation 
for Formation Flying of Spacecraft”, AIAA Guidance, Navigation, and Control 
Conference, Denver, CO, Aug. 14-17, 2000. 


[4] 


H. Khalil, Nonlinear Systems, Prentice-Hall, Inc., 1996, pp. 115-116 and 626-629. 




[5] M. Krstic, I. Kanellakopoulos, and P. Kokotovic, Nonlinear and Adaptive Control 
Design, John Wiley & Sons, Inc., 1995, p. 491. 

[6] O. Egland and J.M. Godhaven, “Passivity-Based Adaptive Attitude Control of a 
Rigid Spacecraft”, IEEE Transactions on Automatic Control, Vol. 39, No. 4, 
April 1994, pp. 842-846. 

[7] B. Wie, Space Vehicle Dynamics and Contrl, AIAA Education Series, 1998, pp. 
318-323. 

[8] R.M. Sanner, “Adaptive Attitude Control Using Fixed and Dynamically 
Structured Neural Networks”, AIAA Guidance, Navigation and Control 
Conference, San Diego, CA, July 29-31, 1996. 




